Experimental determination of Ramsey numbers with quantum annealing 
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Ramsey theory is a highly active research area in mathematics that studies the emergence of order 
in large disordered structures. It has found applications in mathematics, theoretical computer sci- 
ence, information theory, and classical error correcting codes. Ramsey numbers mark the threshold 
at which order first appears and are notoriously difficult to calculate due to their explosive rate of 
growth. Recently, a quantum algorithm has been proposed that calculates the two-color Ramsey 
numbers R(m,n). Here we present results of an experimental implementation of this algorithm 
based on quantum annealing and show that it correctly determines the Ramsey numbers -R(3, 3) 
and R(m, 2) for 4 < m < 8. The R(8, 2) computation used 84 qubits of which 28 were computational 
qubits. This computation is the largest experimental implementation of a scientifically meaningful 
quantum algorithm that has been done to date. 

PACS numbers: 03.67.Ac,02.10.Ox,89.75.Hc 



In recent years first steps have been taken towards the 
goal of experimentally realizing the computational ad- 
vantages promised by well-known quantum algorithms. 
As with any nascent effort, these initial steps have 
been limited. To date the largest experimental imple- 
mentations of scientifically meaningful quantum algo- 
rithms have used just a handful of qubits. For circuit- 
based algorithms, seven spin-qubits were used to factor 
15; while for adiabatic quantum algorithms, four spin- 
qubits were used to factor 143. In both cases specialized 
versions of the quantum algorithms were needed to allow 
factoring with such small numbers of qubits. Although 
factoring was the focus of both experiments, other scien- 
tifically significant applications exist. 

In Ref. 3 a quantum algorithm was proposed for de- 
termining the notorious two-color Ramsey numbers. Ex- 
perimental implementation of this algorithm offers sig- 
nificant challenges as it requires, even for the smallest 
Ramsey numbers, a considerable number of non-local, 
high-order qubit couplings. Because of these challenges, 
Ramsey numbers are an excellent testbed for quantum 
hardware. 

Ramsey numbers are part of an active research area 
in mathematics known as Ramsey theory ! whose cen- 
tral theme is the emergence of order in large disordered 
structures. Problems in Ramsey theory often have an 
intuitive character that is easily grasped. For example, 
imagine that we decide to throw a party. We would like 
the party-size N to be large enough that either there will 
be a group of m people who all know each other (and can 
get the party started), or a group of n people who are 
mutual strangers (and so must mingle, again getting the 
party started). Using Ramsey theory ' , it can be shown 
that a threshold value R(m, n) exists so that as long as 
the party-size N is not less than R(m,n), every party of 
N people will contain either m mutual acquaintances or 
n mutual strangers. The threshold value R(m, n) is an 
example of a two-color Ramsey number. Other types of 



Ramsey numbers exist, though we will focus on two-color 
Ramsey numbers here. 

Ramsey numbers grow extremely quickly and are no- 
toriously difficult to calculate. For (two-color) Ram- 
sey numbers R(m, n) with m,n > 3, only nine are 
presently known'. To dramatize the difficulty of calcu- 
lating R(m,ri), Spencer recounts a famous vignette due 
to Paul Erdos' in which aliens have invaded earth and 
threaten to destroy it in one year unless humans can cor- 
rectly determine R(5,5). Erdos noted that this could 
probably be done if all the best minds and most power- 
ful computers were focused on the task. In this case, the 
appropriate response to the alien demand would be to get 
busy. However, Erdos pointed out that if the aliens had 
instead demanded i?(6, 6) (viz. incrementing m and n by 
1), earthlings should immediately launch a first-strike on 
the aliens as there is little hope of meeting their demand. 

In the Ramsey number quantum algorithm (RQA) 
the calculation of R(m, n) is formulated as an optimiza- 
tion problem which is solved using quantum adiabatic 
evolution. ' Here we present results of an experimental 
implementation of the RQA based on quantum anneal- 
ing and show that it correctly determines the Ramsey 
numbers -R(3, 3) and R(m, 2) for 4 < m < 8. The experi- 
mental computation of R(8, 2) used a total of 84 qubits of 
which 28 were computational qubits, and applied an ef- 
fective interaction coupling 28 qubits. To the best of our 
knowledge, this is the largest experimental implementa- 
tion of a scientifically meaningful quantum algorithm. 

We begin this report with a brief description of the 
RQA. We then discuss the experimental details of imple- 
menting the RQA on a chip of 106 superconducting flux 
qubits using quantum annealing, and follow with a pre- 
sentation of our experimental results. Finally, we close 
with a summary of what has been found. 

1. Quantum algorithm: To start, notice that our 
TV-person party problem can be represented by an TV- 
vertex graph. We associate a vertex with each person, 
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(a) N = 5 (b) N = 6 



FIG. 1: Two-colored graphs and the party problem. 

These 5- and 6-vertex graphs represent instances of parties 
with 5 and 6 people. Pairs of acquaintances are joined by 
blue edges, and pairs of strangers by red edges, (a) This 
graph with N = 5 has no red or blue 3-cliques (i.e. triangles) 
corresponding to a party where no group of 3 people are mu- 
tual friends or strangers, (b) This graph with N = 6 contains 
two blue triangles (with yellow interiors) corresponding to two 
separate groups of 3 mutual friends. Ramsey theory guaran- 
tees the impossibility of coloring the edges of such 6-vertex 
graphs without creating at least one red or blue triangle. 

and an edge is drawn between a pair of vertices only when 
the corresponding people know each other. For m mu- 
tual acquaintances, an edge will connect any pair of the 
m corresponding vertices, while for n mutual strangers, 
there will be no edge between any two of the n corre- 
sponding vertices. The m vertices are said to form an 
m-clique and the n vertices an n-independent set. Our 
party problem is now a statement about graphs: every 
graph with TV > i?(m, n) vertices contains an m-clique 
or an n- independent set. Note that any N- vertex graph 
can be converted into a fully-connected TV-vertex graph 
with two-colored edges by replacing each edge with a blue 
edge and each missing edge with a red edge. This trans- 
forms m-cliques into blue m-cliques and n-independent 
sets into red n-cliques (see Fig. 1). 

We briefly describe the construction of the RQA (see 
Ref. 3 and the Supplementary Information (SI) for a de- 
tailed presentation). As computation of R(m,n) is in- 
timately connected with the presence/absence of edges, 
we associate a bit-variable a v y with each pair of vertices 
{v,v') in an A-vertex graph G. Bit-variable a VtV i = 1 
(0) when v and v' are (are not) joined by an edge in 
G. There are thus L N = (£) = N(N - l)/2 bit- 
variables which we collect into the bit-vector (bit-string) 
a = (02,1 • • • , ajv,i, 03,2, • • • , auv,2, • • • , ajv,jy-i) of length 
Ln. Notice that an A- vertex graph G determines a 
unique bit-string a, and vice versa . 

Ref. 3 (and the SI) showed how to count the number 
of m-cliques (a) and n-independent sets I„ (a) in an 
TV-vertex graph G using its associated bit-string a. Thus 
for any (A- vertex) graph G, h^ n (a) = C^(a) + l£ (a) 
is the total number of m-cliques and n-independent sets 
contained in G. It follows that h^ n (fl) must be greater 
than or equal to zero, and h^ n (a) = if (and only if) 
G a has no m-cliques or n-independent sets. 

The next step in construction of the RQA is to use 



hjn n ( a ) as t ne cos t function for the following combina- 
torial optimization problem (COP): For given integers 
(TV, m, n) and cost function h^ n (a) defined as above, find 
an TV-vertex graph G 0k that yields the global minimum 
of n (a). Notice that if TV < R(m,n), the global min- 
imum is /i^ n (a*) = since Ramsey theory guarantees 
that a graph exists that has no m-clique or n-independent 
set. Furthermore, if TV > i?(m,n), Ramsey theory guar- 
antees that h(a+) > 0. If we begin with A < R(m,n) 
and increment A by 1 until we first find „(<**) > 0, 
then the corresponding A will be exactly R(m,n). We 
will return to this observation shortly. 

Adiabatic quantum optimization (AQO) exploits the 
adiabatic evolution of a quantum system to solve a COP, 
and quantum annealing (QA) provides a physical means 
of implementing AQO. The resultant QA algorithm is 
comparable to simulated anncaling s , but relies on quan- 
tum fluctuations instead of thermal fluctuations. The 
QA algorithm uses the COP cost function to define a 
problem Hamiltonian Hp whose ground-state eigenspace 
contains all COP solutions. The QA algorithm evolves 
the state of an Ljv qubit register from the ground-state of 
an initial Hamiltonian Hi to a ground-state of Hp with 
probability approaching 1 in the adiabatic limit. The 
time-dependent Hamiltonian H(t) driving the QA algo- 
rithm is H(t) = A(t/t f )Hi+ B(t/tf)H P , where tf is the 
algorithm run-time, and adiabatic dynamics corresponds 
to tf — > oo. Here A(t/tf) (B(t/tf)) is a positive mono- 
tonically decreasing (increasing) function with A(l) = 
(5(0) = 0). 

To convert the Ramsey COP into a AQO/QA com- 
putation, each of the bit-variables a v y is promoted to 
a qubit, thus associating a qubit with each vertex pair 
(v,v'). The bit-strings a now label the computational 
basis (CB) states \a), and the Ramsey problem Hamilto- 
nian Hp is defined to be diagonal in the CB with eigen- 
value h!^ n (a). By construction, the ground-state energy 
of Hp vanishes if (and only if) there is a graph with no 
m-cliques or n-independent sets. The initial Hamiltonian 
Hi is the standard one for AQO and appears in the SI. 
Its unique ground-state is the uniform superposition of 
all CB states which is easy to prepare experimentally. 

The quantum algorithm for computing i?(m, n) begins 
by choosing A to be a strict lower bound of the Ram- 
sey number: A < R(m,n). The QA algorithm is run 
on Ln qubits and the energy E is measured at the end 
of algorithm execution. In the adiabatic limit the result 
will be E = since A < i?(m,n). The value of A is 
incremented A — > A + 1, the QA algorithm is re-run on 
Ljv+i qubits, and the final energy is measured. This pro- 
cess is repeated until E > first occurs, at which point 
the current value of A will be equal to R(m,n). In any 
real application of the QA algorithm the evolution will 
only be approximately adiabatic. Thus the probability 
the measured energy will be the ground-state energy will 
be 1 — e. To circumvent this, the QA algorithm must be 
run k ~ C(ln[l — 5]/lne) times so that with probability 
8 > 1 — e, at least one of the measurement outcomes will 
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be the ground-state energy. By making k large enough, 
5 can be brought arbitrarily close to 1. 

For experimental implementation of Ramsey AQO the 
bit-vector a must be transformed to the spin-vector s = 
2a — 1 so that a v y = 1 (0) maps to s v>v i — 1 (—1). In a 
slight abuse of notation we write C£ (s) = C^((l+s)/2), 
I?(8) = J^((l + «)/2)), and = C%{s) + 1» (s). 

2. Experimental constraints: Our hardware im- 
plementation of QA uses RF SQUID flux qubits. Each 
qubit is a superconducting loop interrupted by Joseph- 
son junctions, and the states |0) and |1) correspond to the 
two directions of circulating current about the loop." The 
chip hardware uses Josephson-junction-based devices to 
produce pairwise qubit-coupling. 10 The low energy dy- 
namics of the chip is well-described by a quantum Ising 
model in a transverse held with Hamiltonian: 

As the qubit couplers are physical devices with defi- 
nite locations on the chip, there are constraints on the 
qubit coupling/connectivity. Fig. 2(a) shows the layout 
of qubits and couplers on the chip. The local biases {hi} 
and coupling strengths { Jy} may be programmed to val- 
ues in the ranges [—2,2] and [—1, 1] respectively, and the 
experimentally measured functional forms of the inter- 
polation functions A(t/tt) and B(t/tt) appear in the SI. 
For further details of the chip hardware, see Ref. 11. 

The cost function h^ n {s) is still not ready for experi- 
mental implementation for two reasons: (a) there are k- 
spin interactions with k > 2; and (b) the spin couplings 
do not correspond to the qubit couplings in Fig. 2(a). 
These obstacles are removed as follows. 

(a) Reduction to pairwise coupling: The SI shows that 
C„\{s) involves interactions coupling (™) spins, while 
In i s ) couples (2) spins. These interactions must be re- 
duced to pairwise coupling if (j) is to be realized 
experimentally. We illustrate how such a reduction can 
be achieved by reducing the 3-bit coupling term 010203 
to pairwise coupling using an ancillary bit- variable b and 
penalty function 

P(ai, a 2 ; 6) = aya-z — 2(ai + a 2 )b + 3b. (1) 

Notice that P(oi, a 2 ; b) = (> 0) when the input values 
for 01, 02, and b satisfy b — a\a,2 (b 7^ 0102). Now con- 
sider the quadratic cost function h(b) = ba 3 +fiP(aia,2; b) 
for given values of [i and a%. For \x sufficiently large, 
h(b) is minimized when b satisfies the equality constraint 
b = a\02 which causes the penalty function to vanish. 
The optimum cost is then h(b = aia 2 ) = 010203 which 
reproduces the 3-bit coupling term using a quadratic cost 
function. This example is generalized in Sec. 3 to produce 
the quadratic cost function used to calculate R(m, 2). 

(b) Matching spin to qubit connectivity: A cost func- 
tion with only pairwise spin coupling may still not be ex- 
perimentally realizable as the spin couplings needed may 
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FIG. 2: (a) Layout of qubits and couplers. The chip 
architecture is a 4 x 4 array of unit cells, with each unit cell 
containing 8 qubits. Within a unit cell, each of the 4 qubits 
in the left-hand partition (LHP) connects to all 4 qubits in 
the right-hand partition (RHP), and vice versa. A qubit in 
the LHP (RHP) also connects to the corresponding qubit in 
the LHP (RHP) of the units cells above and below (to the 
left and right of) it. Qubits are labeled from 1 to 128, and 
edges between qubits represent couplers with programmable 
coupling strengths. Grey qubits indicate usable qubits, while 
white qubits indicate qubits that could not be calibrated to 
operating tolerances and were not used. All experiments were 
done on a chip with 106 usable qubits. (b) R(8, 2) embed- 
ding for qubit connectivity. Embedding used to compute 
R(8, 2) that produced the needed qubit couplings. In low 
energy states like-colored qubits have the same Bloch vector 
and constitute a single effective qubit. This allows an indirect 
coupling of qubits that are not directly coupled in hardware. 



not match the qubit couplings available on chip. The 
primal graph (PG) of a quadratic cost function is the 
graph whose vertices are the spin variables, and whose 
edges indicate pairwise-coupled spins. An arbitrary PG 
can be embedded into a sufficiently large qubit graph 
having the structure of Fig. 2(a). An embedding maps 
a PG vertex to one or more vertices in the qubit graph, 
where the image vertices form a connected subgraph of 
the qubit graph. We link this connected set of qubits 
together with strong ferromagnetic couplings so that in 
the lowest energy state these qubits have identical Bloch 
vectors. For example, to couple qubits 104 and 75 in 
Fig. 2(a) (which are not directly coupled) with coupling- 
strength J , we ferromagnetically couple qubits 104, 112, 
and 107 using strongly negative Jio4,ii2 and Jio7,ii2 val- 
ues. The desired coupling is then carried along the edge 
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FIG. 3: Primal graph for C™(a,6). The penalty func- 
tion P(a;b) imposes the constraints = a,L~\a.L and 
bj = ajbj+i. The primal graph displays the pairwise cou- 
plings appearing in C™(a;b) = ai&2 + fJ.P(a;b). 



connecting qubits 107 and 75 with 775,107 = 3 ■ 

3. Ramsey energy functions: We considered a 
number of Ramsey problems which could be solved using 
the 106 qubits available in Fig. 2(a). Specifically, R(ra, 2) 
with 4 < m < 8 and i?(3, 3). 

(a) R(m, 2): Since an TV- vertex graph G a with N < m 
cannot contain an m-clique, it follows that C^(a) = 
for all G a . Thus, for N < m, h^ 2 (a) = I? (a) = 
5i + ••• + a,L N , where Hi = 1 — aj. This produces 
a problem Hamiltonian Hp with Ln uncoupled qubits 
which is easily mapped onto the chip. Now consider 
N = m. Defining L — L m — (™), we have C™(a) = 
aia 2 ■ ■ ■ a L ^a L , and /i™ 2 (a) = C™(a) + I? {a). The L- 
bit interaction in C™(a) is reduced to pairwise coupling 
by introducing: (i) ancillary bit variables &2>'" - >&l-1j 
and (ii) imposing the constraints = a^-ia^ and 

6j = a,jbj + i (j = 2, • • • , L — 2) through the penalty func- 
tion P(o;6) = P(az,_i,ai;6i_i) + X}j=2 2 -ffe) & j) 
(see Eq. (1)). The R(m,2) cost function for N = m 
is then /i™ 2 (a,&) = {aib 2 + [J. P(a;b)} + Hp(a), where 
ji = 2 is the penalty weight value used in all R(m,2) 
experiments. Making the substitutions 2o = s a + 1 and 
2b = St, + 1 expresses the cost function in terms of Ising 
spin variables. 

The PG for pairwise interactions present in ft,™ 2 (a,6) 
appears in Fig. 3. We have embedded this PG into the 
hardware up to N — m = 8. In Fig. 2(b) we display 
the embedding that was used to determine R(8,2), and 
which used 28 computational qubits, 26 ancilla qubits to 
reduce interactions to pairwise, and 30 qubits to match 
the PG connectivity in Fig. 3 to the qubit connectivity 
in Fig. 2(b) for a total of 84 qubits. 

(6) i?(3,3): We also determined i?(3,3) by examining 
N = 4,5, 6. The cost functions for these cases are 

^3,3( a ) — /l,2,4 + /l,3,5 + 72, 3,6 + 74,5,6) 

^3,3( a ) — /l,2,S + /l,3,6 + flAJ + /2,3,8 + f 2,4,9 + 

/3,4,10 + /5,6,8 + /5,7,9 + fe,7,10 + /8,9,10) 
^3,3( a ) — A, 2, 6 + /l,3,7 + /l,4,8 + /l,5,9 + /2,3,10 + /2,4,11 
+ /2,5,12 + /3,4,13 + /3,5,14 + /4,5,15 + fa,7,W 
+ h,8,Xl + h,9,12 + /7,8,13 + /7,9,14 + /8,9,15 
+ 710,11, 13 + 710,12,14 + 711,12,15 + 713,14,15) 



where fi.j,k = a-idjCik + aiajCik- Notice that fij t k can be 
re-written as fij.k = —2+ai+aj+ak+a,iaj+aiak+aja>k, 
which only contains pairwise couplings, making ancillary 
6-spins unnecessary. The largest of these problems is 
N = 6 whose PG has 15 vertices and 60 edges. We 
can reduce its size slightly by exploiting the identity 
hm,n( a ) = hnrnfi)- ^ or m = n tms yields a two-fold 
symmetry: if a* is a global minimum of m , so is a*. 
Thus, we can fix one variable (say a\ = 0) and optimize 
over the remaining variables a' . Optimal solutions then 
have the form (0,a^) and (l,a^). With this simplifica- 
tion, the PG of h® 3 (0,a') has 14 vertices and 52 edges. 
The SI shows its embedding into the chip hardware. 

4. Results: To solve a given Ramsey COP specified 
by the parameters h and J, the chip must first be pro- 
grammed to fix these values in hardware. For the largest 
problem we solved (R(8, 2) using 84 qubits) this took 
roughly 270 ms. After programming we iterate many cy- 
cles of annealing and readout. Each annealing cycle has 
duration tf — 1 ms, and readout of the qubits takes 1.5 
ms per sample. Programming only occurs once so the 
total runtime required to obtain S Ramsey output sam- 
ples is 270 + (1 + 1.5)5 ms. As the hardware is an analog 
device, there is limited precision to which h and J can 
be specified. For a COP whose ground state is sensitive 
to parameter settings this could pose serious difficulties. 
However, the Ramsey COP requires specification of only 
a few distinct integral values, and ground states are quite 
stable to parameter perturbations. 

Although most discussions of QA assume zero- 
temperature, all experiments here were performed at 20 
mK, or roughly 0.4 GHz. For comparison, peak values 
of the interpolation functions A(t/tf) and B(tjtf) are of 
order 10 GHz (see the SI). Refs. 12,13 show that finite 
temperature need not destroy the efficacy of QA. In our 
experiments we select the lowest energy configuration ob- 
served over many annealing cycles to compensate for the 
stochastic influence of non-zero temperature. 

Figs. 4 and 5 present our results for R(8, 2) and R(2>, 3). 
The SI contains our results for R(m, 2) with m = 4, • • • ,7 
which could not be included here due to space limita- 
tions. Both Figures display histograms that plot the 
relative-frequency of energy values obtained by program- 
ming the chip and running 10 5 annealing and readout 
cycles, yielding 10 5 s-spin configurations. In the main 
figures, histograms of the energies of the Ising problem 
sent to the hardware are plotted. These Ising cost func- 
tions include the Ramsey cost function n (s), and the 
ferromagnetic penalties enforcing the equality constraints 
amongst qubits that represent the same PG spin vari- 
able. The ferromagnetic penalty weight A was adjusted 
so that at least 85% of output s-configurations satisfied 
the equality constraints. These feasible spin configura- 
tions were translated back to the original a variables and 
the cost /energy function h^ n (a) evaluated. The result- 
ing energy values were binned and plotted in the inset 
histograms. 

Fig. 4 presents our results for i?(8, 2) which was the 
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shows our results for N — 5 and 6, while N — 4 appears 
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FIG. 4: Energy histograms for R(8, 2) on graphs below 
(N — 7) and at (TV = 8) the Ramsey transition. 

largest of the R(m, 2) problems considered. Of the s- 
conhgurations returned by the hardware for TV = 8, 
approximately 65% are global minima of /i§ 2 - Exam- 
ining the inset histograms for N = 7 (8) we see that: 
(i) h m i n = (1); (ii) the probability for h = (1) is ap- 
proximately 1.0 (0.65); and (iii) the number of optimal a- 
conhgurations/graphs is 1 (29). The energies h min found 
for TV = 7 and 8 agree with the final ground-state (GS) 
energies found in Ref. , indicating that QA hnds the final 
GS with high-probability. As h min jumps from — > 1 as 
N goes from 7 — > 8, the Ramsey protocol correctly 5 iden- 
tifies R(8, 2) = 8. Finally, Ref. 3 showed that the number 
of optimal graphs for N = 7 (8) is 1 (29) which agrees 
with what was found by QA. For N = 7, the unique op- 
timal a-configuration corresponds to the graph in which 
every pair of vertices is connected by an edge and so has 
no 2-independent sets or 8-clique and so has h\ 2 = 0. 
For N — 8, the 29 optimal a-configurations correspond 
to graphs G a with hg 2 (a) = 1 which are the twenty- 
eight 8-vertex graphs containing a single 2-independent 
set, and the single 8-vertex graph containing an 8-clique. 
We also examined #(3,3) with N = 4,5,6. Fig. 5 



FIG. 5: Energy histograms for R(S, 3) on graphs below 
(TV = 5) and at (N — 6) the Ramsey transition. 



in the SI. Together, they show that the minimum energies 
for N = 4, 5, 6 are 0, 0, 2, respectively, and these occur 
with probabilities of approximately 0.93, 0.94, and 0.91. 
The energies h m i n agree with the final GS energies found 
in Ref. indicating that QA again finds the final GS with 
high-probability. As h m i n jumps from — > 2 as N goes 
from 5 — > 6, the Ramsey protocol correctly 1 identifies 
i?(3, 3) = 6. Finally, Ref. showed that the number of 
optimal graphs for N = 4, 5, 6 is 18, 12, 1760, respectively, 
in excellent agreement with the QA results of 18, 12, 1758. 

5. Summary: We presented results of an experimen- 
tal implementation of a quantum algorithm that com- 
putes two-color Ramsey numbers R(m, n) . We showed 
that the algorithm correctly determined R{m, 2) = m for 
4 < m < 8, and i?(3,3) = 6. A summary of our results, 
along with corresponding theoretical predictions, appear 
in Table I. The R(8, 2) computation used 84 qubits of 
which 28 were computational qubits, and to the best of 
our knowledge is the largest experimental implementa- 
tion of a scientifically meaningful quantum algorithm. 
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TABLE I: Results 1 for Ramsey numbers R(m, 2) = m for 4 < m < 8 and i?(3, 3) = 6. 
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N is the number of graph vertices; E ga and D are the ground- 
state energy and degeneracy, respectively, for the problem Hamil- 
tonian Hp; and for each Ramsey number, the experimental results 
are followed by the theoretical predictions from Ref. ! in parenthe- 
sis. 
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The structure of this Supplementary Information (SI) is as follows. We begin in Section I with a description of how 
the Ramsey number quantum algorithm is constructed. Section II then presents those Ramsey number experimental 
results that could not be included in the manuscript due to space limitations. Section II A presents the results for 
i?(3, 3) with TV = 4, while Section II B presents the results for the Ramsey numbers R(m, 2) with 4 < m < 7. For 
easy reference, Section II A also includes the data for i?(3,3) with TV = 5 and 6 which appears in Section 4 of the 
manuscript. Finally, Section III displays the embedding of the Ramsey energy functions into the chip for i?(3, 3) 
with TV = 6 (Section III A) and i?(8, 2) with TV = 8 (Section MB). These embeddings represent the most complex 
embeddings we encountered in our experimental determination of, respectively, diagonal (R(m, m)) and non-diagonal 
(R(m, n), m ^ n) Ramsey numbers. 



I. RAMSEY NUMBER QUANTUM ALGORITHM 

We briefly describe the Ramsey number quantum adiabatic algorithm (see Ref. I for details). We begin by estab- 
lishing a 1-1 correspondence between the set of TV-vertex graphs and binary strings of length L = TV(TV — l)/2. To 
each TV-vertex graph G there corresponds a unique adjacency matrix A(G) which is an TV x TV symmetric matrix 
with vanishing diagonal matrix elements, and with off-diagonal element a^.,- = 1 (0) when distinct vertices i and j 
are (are not) joined by an edge. It follows that A(G) is determined by its lower triangular part. By concatenating 
column-wise the matrix elements cii.j appearing below the principal diagonal, we can construct a unique binary string 
17(G) of length L for each graph G: 

g{G) = a 2 ,i • • • <Xn,\ 03,2 • ■ ■ ajv,2 ■ • • o,n,n-i- (1) 

Given the string g{G), the following procedure determines the number of m-cliques in G. Choose m vertices, 
S a = {vi, ■ . ■ ,v m }, from the TV vertices of G and form the product C a = Yi^vies ) a ""i,'»k- Note that C a = 1 
when S a forms an m-cliquc; otherwise C a — 0. Now repeat this procedure for all p = C(N,m) ways of choosing m 
vertices from TV vertices, and form the sum C(G) = Y)^^ C a . By construction, C(G) equals the number of m-cliques 
contained in G. A similar procedure determines the number of n-independent sets in G. Briefly, choose n vertices 
T a — {t>i, . . • , v n } from the N vertices in G, and form the product T a — Y]^^v k ^T ) ®Vj,vki where a Vj Vk = 1 — a Vj Vk . 
Note that if I a = 1, then T a forms an n-independent set; otherwise l a = 0. Repeat this for all v = C(N,n) ways of 
choosing n vertices from N vertices, then form the sum 1(G) = ^2 a —i^a- By construction, 1(G) gives the number 
of n-independent sets contained in G. Finally, define 

h(G)=C(G)+l(G). (2) 

It follows from the above discussion that h(G) is the total number of m-cliques and n-independent sets in G. Thus 
h(G) > for all graphs G; and h{G) = if and only if G does not contain an m-clique or n-independent set. 

We use h(G) as the cost function for the following combinatorial optimization problem. For given integers N, m and 
n, and with h(G) defined as above, find an TV- vertex graph G* that yields the global minimum of h(G). Notice that 
if TV < R(m, n), the (global) minimum will be /i(G*) = since Ramsey theory guarantees that a graph exists which 
has no m-clique or n-independent set. On the other hand, if TV > R(m, n), Ramsey theory guarantees /i(G*) > 0. If 
we begin with N < R(m,n) and increment TV by 1 until we first find /i(G*) > 0, then the corresponding TV will be 
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FIG. 1: Interpolation functions A(t/tj) and B(t/tf). The functional forms of the experimentally measured interpolation 
functions A(t/tf) and B(t/tf) are shown. The QA run-time tf can be adjusted for times ranging from 20-20 000 /is. For 
comparison, the temperature at which the experiment is performed is shown as the line at roughly 0.4 GHz. 



exactly R{m, n). We now show how this combinatorial optimization problem can be solved using adiabatic quantum 
evolution, which then becomes the basis for a quantum algorithm to compute R(m,n). 

The adiabatic quantum evolution (AQE) algorithm" exploits the adiabatic dynamics of a quantum system to solve 
combinatorial optimization problems. The AQE algorithm uses the optimization problem cost function to define 
a problem Hamiltonian Hp whose ground-state eigenspace encodes all problem solutions. The algorithm evolves 
the state of an L-qubit register from the ground-state of an initial Hamiltonian Hi to the ground-state of Hp with 
probability approaching 1 in the adiabatic limit. An appropriate measurement at the end of the adiabatic evolution 
yields a solution of the optimization problem almost certainly. The time-dependent Hamiltonian H(t) for local AQE 
is 

H(t) = A(t/t f )Hi + B(t/t f )H P , (3) 

where tt is the algorithm run-time; adiabatic dynamics corresponds to tf —± oo; and A(t/tt) (B(t/tf)) is a positive 
monotonically decreasing (increasing) function with ^4(1) = (-B(O) = 0). The experimentally measured functional 
forms of the interpolation functions A[t/tf) and B(t/tf) are shown in Fig. 1. 

To map the optimization problem associated with computing R(m, n) onto an adiabatic quantum computation, we 
begin with the 1-1 correspondence between A^-vertex graphs G and length L = N(N — l)/2 binary strings g(G). From 
Eq. (1) we see that position along the string is indexed by the graph edges (i, j). We thus identify a qubit with each 
vertex pair (i, j), and will thus need L qubits. Defining the computational basis states (CBS) to be the eigenstates of 
cr° ® • •• ® o"^ -1 , we identify the 2 L graph strings g(G) with the 2 L CBS: g(G) -> \g{G)). The problem Hamiltonian 
Hp is defined to be diagonal in the computational basis with eigenvalue h(G) associated with eigenstate \g(G)): 

H P \g(G))=h{G)\g{G)). (4) 

Note that the ground-state energy of Hp will be zero iff there is a graph with no m-cliques or rt-independent sets. An 
operator expression for Hp appears in Ref. I . The initial Hamiltonian Hi is chosen to be 

L-l 
i=0 

where P and ct* are, respectively, the identity and x-Pauli operator for qubit i. The ground-state of Hi is the easily 
constructed uniform superposition of CBS . 

The quantum algorithm for computing i?(m, n) begins by setting N equal to a strict lower bound for R(m, n) which 
can be found using the probabilistic method or a table of two-color Ramsey numbers . The AQE algorithm is run 
on Ln — N(N — l)/2 qubits, and the energy E is measured at the end of algorithm execution. In the adiabatic limit 
the result will be E — since N < R(m, n). The value of N is incremented N — > N + 1, the AQE algorithm is re-run 
on -Ljv+i qubits, and the energy is measured at the end of algorithm execution. This process is repeated until E > 
first occurs, at which point the current value of N will be equal to the R{m,n). Note that any real application of 



TABLE I: Results 1 for Ramsey numbers R(m, 2) = m for 4 < m < 8 and R(3, 3) = 6. 
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TV is the number of graph vertices; E ga and D are the ground-state energy and degeneracy, respectively, for the problem Hamiltonian 
Hp; and for each Ramsey number, the experimental results are followed by the theoretical predictions from Ref. in parenthesis. 



AQE will only be approximately adiabatic. Thus the probability that the measured energy E will be the ground-state 
energy will be 1 — e. In this case, the algorithm must be run k ~ C(ln[l — <5]/lne) times so that, with probability 
6 > 1 — e, at least one of the measurement outcomes will be the true ground-state energy. We can make 5 arbitrarily 
close to 1 by choosing k sufficiently large. 

II. REMAINING RAMSEY NUMBER RESULTS 

Ref. made three predictions that allow theory and experiment to be compared. The first is that the Ramsey number 
R(m, n) is the value of TV at which the global minimum h m i n of the cost function „(a) first becomes non-zero. The 
other two are, respectively, the value of the final ground-state energy E gs = h m i n = min a n (a) and its degeneracy 
D. For reference, Table I from the manuscript is included below which summarizes all our experimental results and 
the corresponding theoretical predictions. Examination of the Table shows that there is excellent agreement between 
the two. 

In Sees. II A and II B we present our results for i?(3,3) and i?(m, 2). For each Ramsey number experiment, 
histograms of the final energies are presented. These energies are determined as follows. At the end of each quantum 
annealing run, the qubits are measured in the computational basis, and the outcome is the final spin variable- 
assignment s. As explained in Sec. 4 of the manuscript, for those spin configurations satisfying the equality constraints, 
the spin configuration s is translated back to the original binary a variables. The energy E(a) is then calculated from 
the Ramsey energy function h^ n (a) given in Sec. 3 of the manuscript, and its value entered into the appropriate 
histogram bin. For each Ramsey number experiment, energy histograms are presented over a range of TV values. As 
in Sec. 4 of the manuscript, for each TV, two histograms are given. The main histogram shows the energy function in 
terms of spins variable s, where the energy function includes the Ramsey energy and the ferromagnetic interactions 
enforcing equality constraints between spins representing the same problem variable. The value of the ferromagnetic 
coupling strength A used in the quantum annealing runs appears in the upper left corner of this histogram. The 
inset histogram corresponds to the a-configurations, and displays the binned Ramsey energies E(a) = n (a). It is 
obtained from the spin configurations (of the main histogram) by filtering out those spin configurations not satisfying 
the equality constraints as these are infeasible (in the optimization sense) spin configurations. Each inset histogram 
contains the following information: (i) the set of observed Ramsey energies E, (ii) the probability (relative frequency) 
for the energy E, and (iii) the number of optimal a-configurations that yielded the minimum energy E = h min . 

Having made these preliminary remarks we examine the remaining Ramsey number results. 

A. R(3, 3) = 6 

As in the manuscript, two histograms are given for each TV: the main histogram corresponds to spin energies that 
include contributions from the Ramsey energy function and the energy penalty functions that enforce the embedding 
of the primal graph vertices into the chip. As described above, the energy in the inset histograms is the Ramsey 
energy n (a). The energy histograms for i?(3,3) appear in Fig. 2. 

Examining the inset histogram for TV = 4 we see that: (i) h m i n = 0; and (ii) 18 distinct a-configurations have 
zero-energy, corresponding to 18 distinct graphs with no 3-cliques or 3-independent sets. From the main histogram 
we see that the probability to find an optimal a-configuration is approximately 0.93. Note that h m i n = is exactly 
the TV = 4 final ground-state energy E gs (tf) = found in Ref. , indicating that quantum annealing (QA) finds the 
final ground-state with high probability. The number of optimal graphs found agrees with that found numerically in 
Ref. 1. 

A similar examination of the main histogram for TV = 5 gives: (i) h m i n — 0; and (ii) there are 12 optimal 
graphs/a-configurations. From the main histogram we see that the probability to find an optimal a-configuration is 
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FIG. 2: Energy histograms for R(3, 3) with N = 4, 5, 6. 



approximately 0.94. Again, h min = equals the N = 5 ground-state energy E gs (tj) = found in Ref. 1, and so QA 
finds the final ground-state with high probability. The number of optimal graphs found experimentally agrees with 
that found numerically in Ref. 1. 

For N = 6, we see that: (i) h m i n = 2; and (ii) 1758 optimal graphs were observed. From the main histogram we 
see that the probability to find an optimal a-configuration is approximately 0.91. As shown in Ref. 1, E gs (tf) = 2 
for iV = 6, and so h m i n — E gs (tf), and QA again finds the final ground-state with high probability. Since N — 6 is 
the first N value for which hmin = E gs (tf) > 0, the protocol for the Ramsey quantum algorithm correctly 1 identifies 
i?(3,3) = 6. Finally, Ref. 1 showed that the number of optimal graphs for N = 6 is 1760 so that QA found all but 
two of the 1760 optimal graphs. 

Table I summarizes all of the above results. 



B. R(m,2) = m 

Here we present our experimental results for R(m, 2) with 4 < m < 7. The results for m = 8 appear in the 
manuscript. Since the R(m, 2) discussion closely parallels that of i?(3, 3), we give a more abbreviated presentation. In 
the manuscript we showed that for N < m, the cost function 2 (a) is linear in a and consequently the spins are un- 
coupled when this cost function is translated to spin variables s a . There is thus no need to introduce ancillary s^-spins, 



1 r 
0.9- 
0.8- 
0.7 
>.0.6- 
n 0.5- 

O 

4 0.4- 
0.3 
0.2- 
0.1 - 

o 



feasible=100.0%, optimal=100.0% 



1 

I 0.8 

CO 

| 0.6 
to 

1 04 

1 0.2 

o 

O 





-1 



distinct optima = 1 



-2 -1 
Energy (hardware) 

(a) N = 3 



1 

0.9 
0.8 
0.7 
, 0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




n 



A=1.0, feasible=98.2%, optimal=96.8% 

distinct optima = 

1 0.8 

03 
£1 

2 0.6 

Q_ 

ra 

I °- 4 

1 0.2 



-13 -12 -11 

Energy (hardware) 

(b) N = 4 



FIG. 3: Energy histograms for i?(4, 2) with N = 3 and 4. 



nor do we need to include ferromagnetic equality penalties in the Ising model. For this reason all spin configurations 
are feasible for N < m, and no A values are recorded in the main histograms. All results below are included in Table I. 

R(4, 2) : Fig. 3 contains the energy histograms for i?(4, 2). Examining the inset histogram for iV = 3 (4), we see 
that: (i) h m i n = (1); (ii) the number of optimal graphs is 1 (7); and (iii) from the main histogram, the probability 
to find an optimal a-configuration is approximately 1.0 (0.97). The minimum energies h min — and 1 for N = 3 and 

4 agree exactly with the corresponding final ground-state energies E gs (tf) found in Ref. I , indicating that QA finds 
the final ground-state with high probability. As /i min jumps from — > 1 as N goes from 3 — > 4, the Ramsey protocol 
correctly J identifies i?(4, 2) — 4. Finally, Ref. 1 showed that the number of optimal graphs for N = 3 and 4 are, 
respectively, 1 and 7, which agrees with the number of optimals found by QA. For N = 3, the unique optimal graph is 
the fully connected 3-vertex graph which has no 4-cliques or 2-independent sets, while for N = 4, the optimal graphs 
are the ( 2 ) = 6 graphs with a single 2-independent set and no 4-clique, and the unique fully connected 4- vertex graph 
which has one 4-clique and no 2-independent set. 

R(5, 2) : Fig. 4 contains the energy histograms for i?(5, 2). Examining the inset histogram for N — 4 (5), we see 
that: (i) h m i n = (1); (ii) the number of optimal graphs is 1 (11); and (iii) from the main histogram, the probability 
to find an optimal a-configuration is approximately 1.0 (0.76). The minimum energies h m i n — and 1 for N — 4 and 

5 agree exactly with the corresponding final ground-state energies E gs (tf) found in Ref. , indicating that QA finds 
the final ground-state with high probability. As h m i n jumps from —> 1 as N goes from 4 — > 5, the Ramsey protocol 
correctly 1 identifies i?(5, 2) = 5. Finally, Ref. 1 showed that the number of optimal graphs for N — 4 and 5 are, 
respectively, 1 and 11, which agrees with the number of optimals found by QA. For N = 4, the unique optimal graph 
is the fully connected 4-vertex graph which has no 5-cliques or 2-independent sets, while for N = 5, the optimal 
graphs are the ( 2 ) = 10 graphs with a single 2-independent set and no 5-clique, and the unique fully connected 

5- vertex graph which has one 5-clique and no 2-independent set. 

R(6, 2) : Fig. 5 contains the energy histograms for i?(6, 2). Examining the inset histogram for TV = 5 (6), we see 
that: (i) h min = (1); (ii) the number of optimal graphs is 1 (16); and (iii) from the main histogram, the probability 
to find an optimal a-configuration is approximately 1.0 (0.82). The minimum energies h m i n = and 1 for N = 5 and 

6 agree exactly with the corresponding final ground-state energies E gs (tf) found in Ref. , indicating that QA finds 
the final ground-state with high probability. As h m i n jumps from —¥ 1 as N goes from 5 — > 6, the Ramsey protocol 
correctly J identifies i?(6, 2) — 6. Finally, Ref. 1 showed that the number of optimal graphs for N = 5 and 6 are, 
respectively, 1 and 16, which agrees with the number of optimals found by QA. For N — 5, the unique optimal graph 
is the fully connected 5-vertex graph which has no 6-cliques or 2-independent sets, while for N = 6, the optimal 
graphs are the ( 2 ) = 15 graphs with a single 2-independent set and no 6-clique, and the unique fully connected 

6- vertex graph which has one 6-clique and no 2-independent set. 
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FIG. 4: Energy histograms for R(5, 2) with N = 4 and 5. 
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FIG. 5: Energy histograms for R(6, 2) with N = 5 and 6. 



R(7, 2) : Finally, Fig. 6 contains the energy histograms for R(7,2). Examining the inset histogram for N — 6 (7), 
we see that: (i) h m i n — (1); (ii) the number of optimal graphs is 1 (22); and (iii) from the main histogram, the 
probability to find an optimal a-configuration is approximately 1.0 (0.66). The minimum energies h min — and 1 for 
N = 6 and 7 agree with the corresponding final ground-state energies E gs (tf) found in Ref. 1, indicating that QA 
finds the final ground-state with high probability. As h mm jumps from — > 1 as N goes from 6 — > 7, the Ramsey 
protocol correctly 1 identifies i?(7, 2) = 7. Finally, Ref. I showed that the number of optimal graphs for N = 6 and 7 
are, respectively, 1 and 22, which agrees with the number of optimals found by QA. For N = 6, the unique optimal 
graph is the fully connected 6-vertex graph which has no 7-cliques or 2-independent sets, while for N = 7, the optimal 
graphs are the Q = 21 graphs with a single 2-independent set and no 7-clique, and the unique fully connected 
7- vertex graph which has one 7-clique and no 2-independent set. 
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FIG. 6: Energy histograms for R(7, 2) with N = 6 and 7. 



III. EMBEDDING PRIMAL GRAPH INTO CHIP HARDWARE 



Fig. 7 reproduces Fig. 2a of the manuscript which is included here for convenience. Recall that the primal graph of 
an optimization cost function consists of: (i) vertices that represent bit variables whose values are to be optimized; 
and (ii) edges that connect pairs of interacting bits. The primal graph of „(s) is the same as the primal graph of 
„(o) as the former is obtained from the latter by the substitution s = 2a — 1. 



A. R(3,3) 



As explained in Section 3 of the manuscript, the cost function h\ 3 (o) has 14 variables since we fixed a\ = 0. Fig. 8 
shows the embedding of the primal graph for the problem h® 3 (0, 02, 03, • • • , 015). Note that the qubit labeled i in 
this Figure corresponds to the binary variable aj+i. Like-colored qubits are connected with ferromagnetic couplings 
of strength A along the like-colored edges. Black edges are used to represent the coupling strengths between primal 
variables. The embedding is situated within the chip so that the top-left qubit labeled as 8 in Fig. 8 corresponds to 
qubit 41 of Fig. 7. Note that many other embeddings into the chip are possible, and we make no claims that this 
embedding uses the smallest numbers of qubits to represent the required connectivity. 



B. R(8,2) 



The Ramsey cost function hg 2 (o) is significantly more complex than i? 3 3 as it involves a product term that couples 
28 spin variables. As discussed in the manuscript this requires the introduction of 26 ancillary 6-variables to reduce 
this product to a sum of pairwise interactions. The resulting primal graph for hg 2( a i^) is shown in Fig. 9 which is 
a reproduction of Figure 3 in the manuscript. For this case there are a total of 54 primal a and b variables. The 
embedding of the primal graph in Fig. 9 into the chip is shown in Fig. 10 which is a reproduction of Fig. 2b in the 
manuscript. Note that an additional 30 qubits are needed to complete the embedding, bringing the total number of 
qubits used in the computation to 84. Like-colored qubits represent a single primal variable, though note that certain 
colors had to be reused. As before black edges carry the primal variable couplings. The topmost qubit labeled 5 in 
Fig. 10 corresponds to qubit 17 of Fig. 7. 




FIG. 7: Layout of qubits and couplers. The chip architecture is a 4 x 4 array of unit cells, with each unit cell containing 
8 qubits. Within a unit cell, each of the 4 qubits in the left-hand partition (LHP) connects to all 4 qubits in the right-hand 
partition (RHP), and vice versa. A qubit in the LHP (RHP) of a unit cell also connects to the corresponding qubit in the LHP 
(RHP) in the units cells above and below (to the left and right of) it. Most qubits couple to 6 neighbors. Qubits are labeled 
from 1 to 128, and edges between qubits indicate couplers which may take programmable values. Grey qubits indicate usable 
qubits, while white qubits indicate qubits which, due to fabrication defects, could not be calibrated to operating tolerances and 
were not used. All experiments were done on a chip with 106 usable qubits. 



C. R(m, 2) for 4 < m < 7 

The procedure for finding the chip embeddings for the Ramsey cost functions h™ 2 with m = 4, • • • , 7 is similar to 
that presented for /i| 2 ■ The resulting embeddings are much simpler and are not reproduced here. 
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FIG. 10: Embedding of the primal graph of Fig. 9 into hardware. Qubits are numbered so that qubits labeled i with 
1 < i < 28 correspond to primal variables at, and qubits labeled 29, • • • ,54 correspond to the ancillary variables 62, • ■ • , 627- 



